Fluctuation-induced tunneling conduction through Ru02 nanowire contacts 
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A good understanding of the electronic conduction processes through nanocontacts is a crucial step 
for the implementation of functional nanoelectronic devices. We have studied the current-voltage 
(7-V) characteristics of nanocontacts between single metallic Ru02 nanowires (NWs) and contacting 
Au electrodes which were pre-patterned by simple photolithography. Both the temperature behavior 
of contact resistance in the low-bias voltage ohmic regime and the I-V curves in the high-bias voltage 
non-ohmic regime have been investigated. We found that the electronic conduction processes in the 
wide temperature interval 1-300 K can be well described by the fluctuation-induced tunneling (FIT) 
conduction theory. Taken together with our previous work (Lin et ai, Nanotechnology 19, 365201 
(2008)) where the nanocontacts were fabricated by delicate electron-beam lithography, our study 
demonstrates the general validity of the FIT model in characterizing electronic nanocontacts. 



CD 



o 



> 
rn 

o 



X 



I. INTRODUCTION 

Nanometer-scale devices based on different kinds of 
(metallic, semiconducting, magnetic, etc.) nanowires 
(NWs) have widely been demonstrated in vari- 
ous potential applications, such as nanoelectronics, 
nano-electromechanical systems, optoelectronics, and 
spintronics.Q To characterize intrinsic NW properties as 
well as to interrogate the performance of a nanodevice, 
how to successfully attach submicron electrical contacts 
to individual NWs is a crucial step. Until recently, the 
charge conduction properties of electronic nanocontacts 
and their possible effects on device operations have of- 
ten been overlooked, in particular in the case of metallic 
NWs.[l-P| Obviously, poor knowledge about such prop- 
erties could lead to a less quantitative description or 
even a misinterpretation of the experimental data. Q Be- 
sides, ignorance of possible degradative effects due to the 
nanocontacts might cause damages to nanoelectronic cir- 
cuitries. In this regard, a good understanding of the elec- 
tronic conduction processes through a nanocontact is of 
both scientific and technological importance. 

In a recent paper (Ref. |4j), we have demonstrated 
that the temperature behavior of electronic contact re- 
sistances, Rc{T), formed at the interfaces between vari- 
ous metallic (Ru02, Ir02, and Sn-doped In2 03_a;) NWs 
and lithography-patterned submicron Cr/Au electrodes 
can be well described by the fluctuation-induced tun- 
neling (FIT) conduction model proposed by Sheng and 
co-workers, p, 0| In Ref. 0, the NWs were first disposed 
onto a Si02-capped Si wafer and the delicate electron- 
beam lithography and lift-off technique was applied to 
make the contacting electrodes in a two-probe configu- 
ration. Under proper conditions, one can intentionally 
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make the nanocontact resistance Rc much larger than 
both the metallic NW resistance and the contacting elec- 
trode resistance. Therefore, the measured device resis- 
tance R{T) ~ 2Rc{T), where we have assumed that R^ is 
the same for the two nanocontacts in a two-probe con- 
figuration because the two nanocontacts had been fab- 
ricated simultaneously and under the same conditions. 
The Rc{T) curves between 1.5 and 300 K had been mea- 
sured in the limit of zero bias voltage, where the current- 
voltage (I-V) characteristics were linear, i.e. ohmic. We 
explicitly demonstrated that Rc increased with decreas- 
ing temperature, but with a rate much slower than that 
would be due to the thermal-activation-type conduction 
(oc exp(r*/T), where T* is a characteristic temperature), 
especially at liquid-helium temperatures. 

In this work, we aim to further test the validity of the 
FIT conduction model by employing an alternative nan- 
odevice contact configuration. Here we fabricate individ- 
ual NW devices with the NWs being directly placed on 
the surfaces of two pre-patterned Cr/Au electrodes which 
were made by the simple photolithography. In particu- 
lar, the I-V curves in both the low-bias voltage regime 
and the high-bias voltage (i.e. non-ohmic) regime have 
been measured. We found that our overall experimental 
results can still be well described by the FIT conduction 
theory, indicating the general importance of the thermal 
fluctuation effect in electronic nanocontacts. [gI,!?!] Our ob- 
servations are reported below. 



II. EXPERIMENTAL METHOD 

The Ru02 NWs used in this study were synthesized b 
the thermal evaporation method, as described in Ref. 
The as-synthesized NWs were dispersed by ultrasonic 
agitation into ethanol solution. Several droplets of the 
solution were dropped onto a Si substrate capped with 
a 500 nm Si02 layer. On the Si02 layer, an array 
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FIG. 1. (a) Current-voltage curves at 1.6, 15, 30, 45, 60, 90, 
and 300 K for device A. The inset shows an expanded plot for 
the I-V curve at 1.6 K. (b) A schematic of a contacting Au 
electrode pair bridged with a single NW. The NW is directly 
placed on the surfaces of the Au electrodes, (c) An SEM 
image of device A. The scale bar is 2 /xm. 

of micron-sized electrode pairs consisting of two Cr/Au 
(w 20/100 nm) fingers separated by 5 /im were pre- 
patterned by the standard photolithography and lift-off 
technique. The substrates were inspected by a scanning 
electron microscope (SEM), and those electrode pairs 
which were bridged with only one Ru02 NW were se- 
lected for I-V curve measurements. Electrical leads from 
the Cr/Au electrode pairs to macro Cu wires were con- 
nected with Ag paste. A Keithley model K-220 current 
source and a model K-617 electrometer were used. For 
measurements from 300 K down to 1.5 K, a standard 
^He cryostat equipped with a calibrated silicon diode 
thermometer was utilized. Figures [Ifb) and [Ijc) show 
a schematic diagram and an SEM image of an individual 
Ru02 NW device, respectively. Two devices are stud- 
ied in this work. Their sample parameters are listed in 
Table H 

We mention in passing that, to acquire a proper Rc 
which is much larger than both the metallic NW resis- 
tance and the contacting electrode resistance, we simply 
cast the NW on top of the pre-fabricated contacting elec- 
trodes. By this method, high-resistance contacts are eas- 
ily obtained, d For semiconducting NWs,[l3 El the 
tunnel barrier is an insulating layer inherently formed 
on the surface of the as-synthesized NW. However, for 
the metallic Ru02 NWs,[l3| the tunnel barrier separat- 
ing the metallic NW and the contacting electrode is ex- 
ternally introduced, [l^ and its exact origin needs further 
investigations. 

III. RESULTS AND DISCUSSION 

Figurellja) shows the measured I-V characteristics for 
the device A at several temperatures between 1.6 and 300 
K, as indicated in the caption to Fig. [TJ It is seen that 
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FIG. 2. Zero-bias resistance as a function of temperature for 
two NW devices, as indicated. The symbols are the exper- 
imental data. The solid curves are the least-squares fits to 
Eq. 11]). The dashed curves are the theoretical predictions of 
the same equation but are plotted by directly substituting the 
7i and Tb values extracted from a{T) (Eq. ([5])). For clarity, 
the data for the device B have been shifted up by multiplying 
a factor of 2. 

the I-V curves are symmetric about zero bias voltage at 
all measurement temperatures. At low bias voltages, the 
I-V curves are linear, i.e. ohmic. However, when the bias 
voltage increases up to a certain level, a nonlinear behav- 
ior appears. As expected, the nonlinearity becomes pro- 
gressively pronounced at lower temperatures. The inset 
of Fig. [IK a) indicates that the ohmic regime persists up 
to about 4 mV at 1.6 K. It should be noted that, even at 
our highest bias voltages and at the lowest measurement 
temperatures, the power dissipation in the NW device 
is still small (<2xl0^^ W), and thus the electron over- 
heating effect can be safely ignored. (The heat mainly 
dissipated at the two nanocontacts.) 

From the I-V curves measured at different tempera- 
tures, the temperature behavior of device resistance R 
can be obtained. The R is defined as the reciprocal of 
the slope of the I-V curve around zero bias voltage. Fig- 
ure [2] shows a double logarithmic plot of the measured R 
as a function of temperature for the devices A and B, as 
indicated. Similar to our previous finding in Ref. 0, one 
sees that the R increases with decreasing temperature 
and eventually saturates at a few degrees of kelvin. [l^ 
As has been demonstrated in Ref. 0, the sum of the indi- 
vidual Ru02 NW resistance and the two Cr/Au electrode 
resistances is ^1.5 Ml, which is ^ 3 orders of magnitude 
smaller than R. Therefore, we may write R{T) ~ 2Rc{T) 
for all temperatures, as mentioned. 

The FIT model addresses the electronic conduction 
properties of a mesoscopic tunnel junction. ^ Q More 
precisely, the FIT theory considers how the electron tun- 
neling process is affected by temperature when the size 
of a tunnel junction formed between two large metallic 
regions becomes sufficiently small. The tunnel junction 



3 



TABLE I. Values of relevant parameters for two Ru02 NW devices. Device A (B) has a NW diameter of ~ 70 (~ 100) nm. 
The junction (nanocontact) area A was estimated from the SEM image and given by the product of the diameter of the NW 
and the average extent of the two Au contacting electrodes lying beneath the NW. The specific contact resistivity is defined 
by pc = RcA and pcO = Pc(r — K). The intrinsic NW resistance for device A (B) is ^--^2 ('-^ 1) kO, and weakly dependent on 
temperature, as determined by the four-probe method in Ref. 1^. The total Au electrode resistances are < 10 Q. 



Device 


Method 


R{300 K) 


Roo 


Ti 


To 


A 


w 


</>o 


Ti/To 


PcO 






(kfi) 


(kD.) 


(K) 


(K) 




(nm) 


(meV) 




{kn fim^) 


A 


RiT) 


37 


17.3 


159 


32.2 


0.166 


15.0 


1.7 


4.94 


401 








16.5 


206 


43.9 




13.7 


1.8 


4.69 


299 


B 


RiT) 


118 


9.9 


1313 


234 


0.138 


10.5 


4.4 


5.61 


374 




/h(Vh) 




20.4 


927 


198 




9.7 


3.6 


4.68 


304 



can be conceptually modeled as a potential barrier with 
a barrier height (j)o, a barrier width w, and a lateral area 
A which is the geometrical area at the closest contact be- 
tween the two large metallic regions. In addition to the 
externally applied bias voltage, 14, the FIT theory takes 
into account the effect of the thermal voltage noise, Vth, 
on the tunneling current across the junction. Sheng and 
coworkers realized that the Vth ~ \/A:bT/C could be large 
due to a small A, where fee is the Boltzmann constant, 
and C (X A is the capacitance of the junction. That is, a 
small junction area A leads to a notable T4h- This large 
Vth is of random sign with regard to V^ and is shown to 
significantly modify the effective tunnel barrier seen by 
the electrons. As a result, the charge transport through 
the junction is profoundly affected. [15| Historically, the 
FIT theory has been formulated to explain the tempera- 
ture dependence of resistance in a wide variety of metal- 
insulator composites and polymers, where a huge number 
of tunnel junctions is involved and an effective junction 
representing the whole network had to be assumed in the 
model. [H-Ul 

According to the FIT theory, [1, 0] at small applied 
electric fields where the I-V characteristic is linear, the 
temperature dependent zero-bias resistance across a sin- 
gle small junction is given by 



R{T)^R^cxp(^-^^^ , 



(1) 



where Rao is a parameter which depends only weakly on 
temperature, and Ti and Tq are characteristic tempera- 
tures defined as 
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and 



To = 
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(3) 



where is the permittivity of vacuum, 27r?i is the 
Planck's constant, and m is the electronic mass. In 
Eq. (HI), k^Ti can be regarded as a measure of the energy 
required for an electron to cross the barrier, and Tq is the 



temperature well below which thermal fiuctuation effects 
become insignificant. Notice that, in the present work, 
the junction (nanocontact) area A is experimentally de- 
termined by SEM and is not an adjustable parameter. 
For comparison, in previous studies of metal-insulator 
composites, polymers, and nanofibers. [l6l - [20j the num- 
ber of tunnel junctions and their individual areas are un- 
known and can not be experimentally determined. 

For large applied electric fields where the I-V curve 
becomes nonlinear, the FIT theory[^ 0] also gives an 
expression for the I-V characteristic of the junction at 
the temperature T: 



/h(Vh,r) = /h,sexp 



-a(T) 1 



Vh_ 

Vh. 



\Vh\< Vh, 



(4) 

where the subscript 'h' denotes the property at 'high' bias 
voltages. The saturation current, Jh,s) and the critical 
voltage, Vh,c, are parameters which depend only weakly 
on temperature. a(T) is a parameter which contains the 
dominant temperature effect on the I-V curves and has 
the form 



a{T) 



To + T 



(5) 



It should be noted that the ohmic regime corresponding 
to small applied voltages is not included in the deriva- 
tion of Eq. (Ill) . Hence, the temperature dependent zero- 
bias resistance (Eq. (|T|)) cannot be directly derived from 
Eq. (HI) by taking the Vh^O hmit of (d/h/dVh)"^ 

Apparently, by measuring the I-V curves at various 
temperatures, one has two different routes to extract the 
values of Ti and Tq. The first route involves the low- 
bias voltage ohmic I-V regime. The two values can be 
obtained by fitting the temperature behavior of R to 
Eq. (|T|). On the other hand, the second route involves 
the high-bias voltage non-ohmic I-V regime. The value 
of a(T) can first be extracted at each temperature by fit- 
ting the measured Ih{Vh) curve to Eq. (|4|) (with the ohmic 
regime being excluded). Then, the values of Ti and To 
can be determined by fitting a(T) to Eq. (|S|). These 
two complementary routes can provide a self-consistency 
check of the phenomenological FIT theory. 

Return to Fig. [5| The symbols are the experimental 
data and the solid curves are the least-squares fits to 
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FIG. 3. Current-voltage characteristics at high bias voltages 
for (a) device A and (b) device B at several temperatures, as 
indicated. The symbols are the experimental data. The solid 
curves are the least-squares fits to Eq. Q. 
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FIG. 4. The parameter a in Eq. ([4]) as a function of tem- 
perature for devices A and B, as indicated. The solid curves 
are the least-squares fits to Eq. ((5]). The dashed curves are 
the theoretical predictions of the same equation but plotted 
by directly substituting the T\ and To values extracted from 
the R{T) fits to Eq. Q. 



Eq. ^ with i?oo, T\ and Tq as adjustable parameters. It 
is seen that the predictions of Eq. ([T]) can well describe 
the experimental data. Our fitted i?oo, Ti, and Tq values 
are listed in Table HI (The dashed curves also describe 
the experimental data reasonably well, see further dis- 
cussion below.) We stress again that the behavior 
shown in Fig. [2] is distinct from the thermal-activation- 
type conduction; it does not obey the K oc exp(T*/r) law 
especially at liquid-helium temperatures (see, for exam- 
ple, the Fig. 1 of Ref. H). The i?cx=, Ti, and Tq values 
are in accord with those corresponding values obtained 
in Ref. % 

Figures [Slja) and [3jb) show representative I-V curves 
for the two devices in the nonlinear regime and at several 
temperatures, as indicated. The symbols are the experi- 
mental data and the solid curves are the least-squares fits 
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FIG. 5. Variation of saturation current /h,s with temperature 
for (a) devices A and (b) device B, and variation of critical 
voltage Vh,c with temperature for (c) devices A and (d) de- 
vice B. Notice that these two parameters are temperature 
independent within our experimental uncertainties. 



to Eq. (HI) with /h^s, Vh^c, and a(T) as adjustable parame- 
ters. Clearly, the predictions of Eq. Q can well describe 
the experimental data in this non-ohmic regime. Our fit- 
ted values of a as a function of temperature is plotted in 
Fig. m and our fitted values of /h,s and l^_c as a function 
of temperature are plotted in Fig. [5l Furthermore, the 
extracted values of a(T) have been least-squares fitted to 
Eq. (O, as indicated by the solid curves in Fig. The 
Ti and Tq values thus inferred are listed in Table HI 

We now examine the self-consistency of the FIT the- 
ory by comparing the T\ and Tq values extracted from 
the two different manners, i.e. i?(r) and /h(Vh). To do 
this, we plot in Fig. [2] the predictions of Eq. ([T|) by di- 
rectly substituting in the T\ and Tq values inferred from 
a(T) while leaving R^o as the sole adjustable parameter. 
The results are indicated by the dashed curves in Fig. [2] 
Clearly, the experimental data are fairly well reproduced 
in this manner. The i?oo value thus extracted is listed in 
Table [J On the other hand, we now substitute those T\ 
and To values extracted from R(T) (Eq. ^) into Eq. © 
to calculate a(r). The calculated results are plotted as 
the dashed curves in Fig. 01 which is seen to well describe 
the a(T) value for device A. (For device B, the discrep- 
ancy is somewhat larger, due to the lack of experimental 
data points at temperatures below 100 K.) In short, our 
T\ and Tq values inferred from the two methods differ by 
an amount of ~30%. This amount of difference is satis- 
factorily acceptable, and is somewhat smaller than that 
(~ 40%) reported in a previous study of carbon nanotube 
network. |2ll| 

In the FIT theory, the barrier shape in the presence of 
an electric voltage across the junction is modeled by the 
image-force corrected rectangular barrier potential origi- 
nally calculated by Simmons. [2^ It was found that differ- 
ent amounts of image-force correction result in different 
shapes of the barrier. With an increase in the electric 
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voltage V {V = Va + Vth is the total voltage), the barrier 
is effectively lowered and narrowed. Hence, there exists 
a (temperature independent) critical voltage, Vh,c, above 
which the peak value of the barrier is fully suppressed. 
Since the barrier is no longer visible to the electrons when 
V > Vh^c, a further increase in V would only increase the 
energy of the electron but leave the barrier transmission 
factor unchanged (the transmission factor now becomes 
unity). This leads to a saturated tunneling current, /h,s, 
which is insensitive to the voltage ^ for V > Vh^- Be- 
sides, it has been shown by numerical evaluation^ that 
/h,s depends only weakly on temperature. Therefore, the 
I-V curves at all measurement temperatures are expected 
to saturate to a constant /h,s- Indeed, a signature of such 
saturation is evident in Fig. [3] in both devices. 

Figure [S] shows the variation of our fitted /h,s and 
14i,c with temperature for the devices A and B. It is 
seen that both parameters are essentially temperature 
independent, as expected by the FIT theory. We may 
evaluate the critical electric field, i?h,c, at the junction 
(nanocontact) from T4i,c for the device A. Recall that 
the voltage mainly dropped at the two nanocontacts (be- 
cause i?~2i?c)- Taking a value Vh,c«43 mV (Fig. El 
and dividing it by 2w w 30 nm (Table |T| , we obtain 
i;h,cWl.4xlOS V/m. Theoretically, the FIT model pre- 
dicts a critical electric field given by £'h,c = 4(/'o/e'i«.@ 
(This corresponds to a critical voltage a few times of 
the junction barrier height.) Substituting the 0o and 
w values from Table U we obtain iJh,c~5xlO^ V/m. 
This agreement again supports the self-consistency of 



the FIT model. We mention in passing that the mi- 
croscopic processes of the FIT conduction mechanism 
have recently been theoretically reexamined by consid- 
ering the electronic wave transmission through long, nar- 
row channels. [13 

IV. CONCLUSION 

In two-probe individual metallic NW devices and un- 
der appropriate conditions, the measured device resis- 
tance is dominated by the two nanocontact resistances. 
These nanocontact resistances are found to be well de- 
scribed by the fiuctuation-induced tunneling conduction 
theory. The parameters characterizing the nanocon- 
tacts can be quantitatively extracted. Taken together 
with our previous and complementary work (Ref. Ij) , our 
study demonstrates that the physical properties of elec- 
tronic nanocontacts can be reliably inferred from the FIT 
model. This information is useful for the characteriza- 
tions and implementation of nanoelectronic devices. 
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